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ABSTRACT 

We have tested the effect of spatial gradients in stellar mass-to-light ratio (T) on measurements 
of black hole masses (M.) derived from stellar orbit superposition models. Such models construct a 
static gravitational potential for a galaxy and its central black hole, but typically assume spatially 
uniform T. We have modeled three giant elliptical galaxies with gradients a = dlog(T)/(ilog(r) 
from —0.2 to +0.1. Color and line strength gradients suggest mildly negative a in these galaxies. 
Introducing a negative (positive) gradient in T increases (decreases) the enclosed stellar mass near 
the center of the galaxy and leads to systematically smaller (larger) M, measurements. For models 
with a = —0.2, the best-fit values of M. are 28%, 27%, and 17% lower than the constant-T case, 
in NGC 3842, NGC 6086, and NGC 7768, respectively. For a = +0.1, M. are 14%, 22%, and 17% 
higher than the constant-T case for the three respective galaxies. For NGC 3842 and NGC 6086, this 
bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative 
(positive) gradients in T cause the total stellar mass to decrease (increase) and the dark matter 
fraction within one effective radius to increase (decrease). 



1. INTRODUCTION 

Dynamical measurements of black hole masses (M») at 
the centers of nearby galaxies have opened new avenues 
for studying galaxy evolution, by exposing correlations 
between M, and various host galaxy properties. Ongo- 
ing studies of these black hole scaling relations explore 
whether the growth and quenching of black holes and 
galaxies are causally linked, or only loosely connected by 
broad trends in cosmic evolution. One way to improve 
on the determination of the black hole scaling relations 
and the cosmic scatter in M, is to reduce uncertainties 
and systematic biases in measurements of M, in local 
galaxies. 

Among the approximately 70 nearby galaxies with 
dynamically determined M m , about 50 M, have been 
measured using s tellar kinematic data an d orbit super- 
position models (McConnell & Ma 20131 and references 
therein). A number of assumptions made in the models 
can contribute to potential biases in the inferred black 
hole masses. Recent advances include modeling stel- 
lar mass and dark matter as separate components (e.g., 
iGebhardt fc Thom as 2009), and accounti ng for triaxial- 
ity (e.g.. Ivan den Bosch fc de Zeeuw|[20Tol) . 

Another assumption that has been used in most pre- 
vious orbit models is that the stellar mass-to-light ratio, 
T, is constant throughout a galaxy. This assumption is 
clearly a simplification since spatial gradients have been 
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observed for many galaxy properties: e.g., color, metal- 
licity, age, and a-element enhancement, which can all be 
associated with spatial variations in T. Early- type galax- 
ies are often bluer towards their outer regions. Decreas- 
ing metallicity outward is likely the domin ant cause for 
the color gradients in the se galaxies (e.g., iStrom et al.l 
[1971 iTamuraet al|[2000l ). with the age gradient play- 
ing a minor or insignifi cant role (e.g.. iSaglia et al.ll2000t 
iTamura fc Ohtal 120041 ) . This picture is supported by 
recent photometric studies of lar ge statistical s amples 
of early-type galaxies in SDSS (jTortora et al.l l2010a|) 
as well as spectroscopic studies of individual galaxies 
dSpolaor et al.ll2009HRawle et al.ll2010l:lKuntschner et al.l 
120101 ). The spatial coverage of these studies is typi- 
cally limited to the central regions of the galaxies within 
one effective radius (r c ff). Measurements of ten early- 
type galaxies have been extended to larger radii us- 
ing deep observations with integral field spectrographs 
(|Weiimans et al.ll2009t fGreene et al.ll2012ft . The central 
gradients in metallicity and absorption line strengths 
are seen to continue to 2-4 r c g. Gradients in color 
and inferred metallicity in galaxies' globular cluster sys- 
tems persist to ~ 8 r c g (jArnold et al.H20lTI : iForbes et al.l 

MM- 

Tight positive correlations between T and color have 
been found in both late- type and early-type galaxies 
(e.g..lBell fc de Jongj[200TI : iBell et al.ll2003l : ITortora et all 
120 111 ), with redder colors corresponding to larger T. By 
fitting synthetic spec tral models to the observed op- 
tical color gradients, ITortora et al.l (|2011[) have found 
mild negative gradients in T for local early-type galaxies 
with stellar ages > 6 Gyr, over a large range of stellar 
masses (M* ~ 10 9 -10 11 3 M©). Color and mass pro- 
files have also been recently invest igated in galaxies at 
0.5 < z < 2.5 (jSzomoru et al.ll2013h . 

In this letter, we assess the amount of systematic er- 
ror in M m resulting from the standard assumption of 



constant T in the orbi t models. We modif y the orbit 
superposition code by iGebhardt et al.l (|2000f ) and intro- 
duce a spatial gradient in T to the stellar component 
of the potential. We investigate the effects on M, for 
three giant elliptical galaxies that we have previously an- 
alyzed (assuming constant T): NGC 384 2, NGC 6086, 
and NGC 7768 (IMcConnell et al.1 [2011allbl . 12011) . All 
three are massive ellipticals and brightest cluster galax- 
ies (BCGs) with high stellar mass M* and velocity dis- 
persion cr: M* = 1.55 x 10 12 , 1.43 x 10 12 , 1.16 x 10 12 M 
and cr = 270, 318, 257 km s" 1 for NGC 3842, NGC 6086, 
NGC 7768, respectively. The black hole masses were pre- 
dominantly constrained by stellar k inematics from the 
integral-field s pectrographs OSIRIS dLarkin et al.ll2006fl 
and GMOS (jAllington- Smith et al.l l2002t iHook et al.l 
I2004D Stellar kinematics at l arge radii were r ecorded with 
the Mitchell Spectro graph (iHill et all I2008D or adopted 
from the literature (jCarter et al.l 119991 lLoubser et al.l 
[2008h . 

A few prior studies of M, have explored variations in 
T. Recent investigations of SO galaxies NGC 1332, NGC 
3368, and NGC 3489 mo deled bulge and disk components 
with s eparate T values dNowak et al.l 120 10; Rus h et all 
1201 If) . iNowak et al.l (|2007fl assumed a s eparate T value 
for th e nuclear disk in NGC 4486A. IGebhardt et al.l 
(2000) assessed the V - / color gradients in NGC 3379 
and ap plied the corresponding T gradient to one set of 
models; iCappellari et al.l ((20021 ) applied a similar treat- 
ment to IC 1459. For both galaxies, models with a T 
gradient and models with uniform T yielded statistically 
consistent measurements of M,. Here we extend from 
prior investigations by examining systematic trends in 
measured M. over a range of T gradients. 

2. STELLAR ORBIT MODELS AND T GRADIENTS 

We use the axisymmetric o rbit modeling algorit hm of 
Gebhardt etldl (|2000| . l200l . I Thomas et al.l ((20051) , and 
Siopis et all ((20091 ) . In this model, a galaxy is described 
by the density profile 



p{r, 6) = T(r) i/(r, 6) + M.S(r) + Pdm (r) 



(1) 



where v(r, 9) is the luminosity density derived from 
the galaxy's deprojected surface brightness profile, and 
Pdm(7") represents a spherical NFW or cored logarithmic 
(LOG) dark matter profile. Equation ([1]) is used to con- 
struct a gravitational potential, and kinematic models 
of the galaxy are constructed by propagating test par- 
ticles through the potential and computing their time- 
averaged velocities ("orbits") throughout a polar grid. 
Orbital weights are varied to fit the observed kinematics 
of the galaxy but are constrained such that the sum of 
weighted orbits exactly reproduces v(r,6), which is de- 
coupled from variations in T(r). The model velocity dis- 
tributions are then compared to the observed kinematics, 
yielding a goodness-of-fit statistic \ 2 ■ Each instance of 
the model adopts a single density profile and outputs a 
single x 2 value. The best-fit values and confidence limits 
for M. and other parameters in Equation ([IJ are deter- 
mined by analyzing the distribution of \ 2 from many 
models. 

Previous investigations using orbit superposition mod- 
els typically have set T(r) to a uniform value, Tq. To 
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Fig. 1. — Color vs. radius in three giant elliptical galaxies based 
on photometric data from the Hubble Space Telescope and Kitt 
Peak National Observatory. In each panel, the black dashed line 
represents the best log-linear fit to all data points, and the solid 
grey line represents the best log-linear fit to data with r < r e g. For 
NG C 6086, our V-I data extend to only 0.5 r cff (17")- B-iJdata 
from Pcleticr et al. (1990) extend to 1.2r e ff and have a logarithmic 
slope of —0.06. 



(2) 



assess the impact of a gradient in T, we set 

_ dlogT 
dlog r 

For each value of a, we sample a two-dimensional grid of 
values for M m and the normalization of T. 

The galaxies' color gradients indicate likely values of 
a. Figure Q] shows the color profile for each galaxy in 
our sample. The color gradients S(B-R)/ 5 (\og 10 r) and 
(5(y- /)/(5(login r) range f rom -0.02 to —0.06. In Figure 
2 of iTortora et al.l (|2011l ), the matching gradients in T 
are « —0.1 to for SDSS galaxies. Their Figure 3 shows 
the T gradients as a function of stellar mass and velocity 



2 



dispersion and indicates different mean vaiues of a for 
galaxies with different stellar ages. Early-type galaxies 
with central age older than 6 Gyr have a ~ —0.2 to —0.1 
over a large range of a, but show a slight upturn toward 
a = at > 1O U M . 

For NGC 3842, we have an additional two hours 
of on-source data from May 2012 observations with 
the Mitchell Spectrograph, allowing for a direct mea- 
su rement of the T gr adient. As described in detail 
in IGreene et al.l (|2012f ). we construct elliptically aver- 
aged radial bins to maximize the S/N in the spec- 
tra. We have a median S/N of 40 per pixel in the 
largest bin, at r w l.lr e ff- Fi rst, weak emission fea- 
tures are iteratively corrected ()Graves fc Faberl [2010^ 
and velocit y dispersions are measu red at each bin us- 
ing pPXF (jCappellari fe Emsellenl I2004D . Then the 
stellar population properties are measured from th e 
Lick_EW+EZ_Ages algorithm ([Graves fc Schiavonl I2008D . 
We use the a— enhanced models of iSchiavonl (|2007fl ~s- 
suming a Salpcter IMF, to determine the luminosity- 
weighted T at each radius. We find the resulting T gra- 
dient to be a = — 0.13±0.30, broadly consistent with the 
color gradient in NGC 3842. Uncertainties in the inner 
bins are dominated by systematic uncertainties in emis- 
sion line and velocity dispersion measurements, while the 
outer bins include additional uncertainties from sky sub- 
traction. 

Gradients in the stellar initial mass function (IMF) are 
potentially an additional source of T gradients within our 
galaxies. The IMF varies with the [Mg/Fe] line ratio and 
the corresponding inferred tim escale for star formation 
(|Conrov fc van Dokkumir2012[ ). Our line index measure- 
ments from Mitchell Spectrograph data of NGC 3842 are 
consistent with no radial gradient in [Mg/Fe]. Therefore 
we do not expect IMF gradients to alter the likely range 
of a. 

Given the uncertainties discussed above, we have cho- 
sen to test a ranging from —0.2 to 0.1. We employ a 
as a spherical, three-dimensional gradient in the models. 
Integrating our model stellar mass profiles along the line 
of sight, we find the projected gradients to be « 80% 
of the three-dimensional a. The corresponding range of 
projected gradients, —0.16 to 0.08, is appropriate with 
regards to the inferred gradients from photometry and 
spectroscopy. 

3. RESULTS 

Table Q] and Figures [2] and |3] illustrate our best- fit val- 
ues of M. and T from modeling NGC 3842, NGC 6086, 
and NGC 7768 with different index a for the T gradient. 
The stellar orbit models optimize the normalization of 
T, while the radial dependence is set by Equation @ for 
fixed a. Each galaxy's distance D, effective radius r or f, 
and black hole radius of influence, r- m i = GM,/a 2 , are 
listed in Table [1] Since r; n f is much smaller than r c g 
(with the ratio ranging from 0.006 to 0.03), we list T at 
three radii (r- m {, 0.1 r c g, and 1 r e g) for comparsion. We 
also list the average luminosity- weighted stellar mass-to- 
light ratio (T), equal to the total stellar mass divided by 
the total luminosity. 

Table Q] and Figure [5] show that within r; n f , there is 
a mild degeneracy between stellar mass and black hole 
mass. As a increases from negative toward positive val- 
ues, T(n n f) decreases, causing the best-fit M. to in- 



TABLE 1 

Results from Stellar Orbit Models 



a M. T(r inf ) T(r cff /10) T(r cff ) <T) 

(1O 9 M ) (T ,b) (T q ,h) (T , h ) (T q ,k) 

NGC 3842 ( r inf = 1.2" ; r eS = 37.8" ; D = 98.4 Mpc) 



-0.2 


7 i +2.8 
'• ± -2.8 


6 4 +0 - 9 


5 1+0-7 


q O+0.4 
° -0.4 


3.5; 


-0.1 


a.7±2:g 




5 1+0-7 
°' i -0.8 


4 1+0-6 


4.3; 


-0.05 


q q+3.3 
M - u -2.5 


5 5+°- 7 


5 2 +0 ' 7 


4 6 +0 ' 6 


4.7; 







5 2+0-7 


5 2+ ' 7 


r. 9+0.7 
°' z -0.7 


5.2" 


0.05 




48 +0.7 


5 1+0-7 


r 7+0.8 
°- ' -0.8 


5.6; 


0.1 


11.2t»;« 


4 5+0-6 


5 0+ ' 7 
a.u_ 7 


6-3±S:i 


6.2" 


0.2 


1 9 4+ 2,2 
lz '^-2.6 


n Q+0.5 
■ s ' a -0.6 


5 +0 ' 6 


7.8t?; 9 


7.5; 



+0.5 
-0.5 
+0.6 
-0.7 
+0.6 
-0.7 
+0.7 
-0.7 
+0.8 
-0.8 
+0.8 
-0.9 
+0.9 
-1.1 



NGC 6086 (r inf = 0.24" ; r cft = 36.8" ; D = 139.1 Mpc) 



-0.2 




6 l+°' 5 


o fi +0.3 


9 9+0.2 
z,z -0.2 


2.8 


-0.1 


3.o±l:i 


r n+0.4 
0.0_q 4 
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4 - u -0.3 


o 9+O.3 
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4 4+0-4 


4 4+0-4 
*' 4 -0.4 


4 4+0-4 


4.4 


0.1 


4.4+i-j 


3 «+02 


4 8+ - 3 
*-°-0.4 


6 0+ ' 4 
D - u -0.5 


5.6 



+0.2 
-0.3 
+0.3 
-0.3 
+0.4 
-0.4 
+0.4 
-0.5 



NGC 7768 (r int = 0.14" ; r eS = 23.1" ; D = 112.8 Mpc) 



-0.2 


i o+°- 5 


C.O_ 1 2 


43+0.7 


9 y+0.4 
z - '-0.4 


2.7 


-0.1 


1 1+0.5 


6 3+ - 8 


4 8+0-6 


n O + 0.5 
J -°-0.6 


3.8 





i 9+O.5 
1 - z -0.3 


5 0+ 06 


5 0+ ' 6 


5 0+ ' 6 


5.0 
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1 4+0-6 


4 0+ - 4 
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°" :> -0.6 


6 7+ ' 7 
°- '-0.8 
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+0.4 
-0.5 
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Best-fit values and errors in M. and T are the median values 
and 68% confidence limits derived from the cumulative likelihood 
method of McConncll ct al. (2011a). Stellar mass-to-light ratios T 
are reported in R band. (T) is the average luminosity-weighted 
value of T over the whole galaxy. 



crease. At large radii the trend with a is reversed, such 
that the best-fit T(r e ff) and (T) both increase as a in- 
creases. The intermediate radius of ~ 0.1r o ff is roughly 
the pivot point at which the best-fit T has similar val- 
ues for different a. Our kinematic data for each galaxy 
fall mainly within 1 r c g and are most thoroughly rep- 
resented at small radii (r ^ 3", or ~ 0.1r e ff), where 
the OSIRIS and GMOS integral-field spectrographs yield 
measurements at multiple position angles. This kine- 
matic coverage drives the normalization of T, such that 
best-fitting models for each value of a have similar en- 
closed masses near 0.1r o ff for a given galaxy. This is 
also the scale on which stars dominate the total enclosed 
mass: much within ~ 0.1r c ff, the black hole begins to 
contribute, whereas the dark matter begins to dominate 
at ~ 0.5 - 2r eff (Figure i]). 

Figure [3] shows a clear positive correlation between a 
and the best-fit M, value. As a decreases from to 
-0.2, the best-fit M. decreases by 28% in NGC 3842 
and 27% in NGC 6086. For these galaxies, the bias in 
Mm over this range in a is comparable to the statistical 
uncertainty in M, for each suite of models with fixed a. 
NGC 7768 exhibits a slightly smaller bias of 17% in M„ 
versus statistical uncertainties ~ 40% from individual 
trials. 

For each galaxy, the value of x 2 f° r the best-fit model 
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Fig. 2. — Contours of x 2 as a function of T(r; n f) and M», for 
orbit superposition models with different slopes for the logarithmic 
T gradient: a = —0.1 (blue dashed), (red solid), and +0.1 (green 
dotted). For each a, the thick and thin contours represent 68% 
and 90% confidence levels (corresponding to A\ 2 = 1.0 and 2.71), 
respectively. 

(for a given a) changes by less than 1.0 when a increases 
from —0.2 to 0. By contrast, \ 2 f° r the best-fit model 
increases by 2.5 as a increases from to +0.2 for NGC 
3842, and \ 2 increases by 1.0 and 2.8 for NGC 6086 
and NGC 7768, respectively, as a increases from to 
+0.1. Models with a positive gradient in T are therefore 
mildly disfavored by our kinematic data. This trend is 
consistent with the range of a inferred from the three 




Fig. 3.— Best-fit M. vs. T gradient a for models of NGC 3842, 
NGC 6086, and NGC 7768. The y-axis is the ratio of the best-fit 
M m value relative to the best-fit value for models with uniform T 
(a = 0). Error bars represent the statistical errors from individual 
trials (Table [TJ, also normalized relative to M,(a = 0). 

galaxies' negative color gradients shown in Figure [TJ 

Figured] shows the relative contributions of stars, dark 
matter, and the black hole to the total enclosed mass as 
a function of deprojected radius in NGC 3842. Similar 
trends are seen for NGC 6086 and NGC 7768. For sim- 
plicity, we have modeled a single dark matter halo for 
each galaxy, for all values of a. As a decreases from to 
—0.2, our kinematic data prefer a model with more (less) 
stellar mass at r <C 0.1 r; n f (r 3> 0.1 r; n f). Since we have 
assumed a fixed dark matter halo, the resulting enclosed 
stellar mass fraction is decreased at large radii. 

If a range of dark matter profiles were sampled, we 
would expect changes in T(r) to alter the mass and/or 
shape of the best-fitting dark matter halo. Introducing 
a negative T gradient would redistribute stellar mass to- 
ward smaller radii, and the models could reconstruct the 
outer mass profile by fitting a more massive halo, ampli- 
fying the trend for dark matter in Figure HJ Our models 
are are not sensitive to the inner slope of the dark matter 
profile, as the enclosed mass at small radii is dominated 
by stars and the black hole. 

Since the dark matter halo and T gradient both cause a 
non-uniform total mass-to-light ratio, one might suspect 
that optimizing the dark matter profile for each value 
of a would decrease the variations in T normalization 
and M.. While our kinematic data for NGC 3842, NGC 
6086, and NGC 7768 do not extend far enough to tightly 
constrain the dark matter profile, we have performed a 
simple test by varying a in models of NGC 3842 with 
no dark matter. Even without dark matter, the overall 
trend of T and M, versus a is similar to those depicted 
in Figures [2h. and [31 

4. DISCUSSION 

We have run stellar orbit superposition models of three 
giant elliptical galaxies, assuming a range of gradients 
defined by a = dlog(T)/dlog(r) in the stellar mass-to- 
light ratio T. We find that the best-fit black hole mass 
M, varies systematically with the strength and sign of 
a. As a decreases from to —0.2, the best-fit values of 
M. decrease by 28%, 27%, and 17% for NGC 3842, NGC 



4 



1.0 



I, J I I I 1 1 1 1 1 — I — I I I I 1 1 1 1 — I I I I 1 1 1 1 1 — I — I I I I I 

' stars J 

•. a = -0 2 black hole . 




Fig. 4. — Enclosed mass fractions vs. deprojected radius for 
stars, dark matter, and a black hole, in the best-fit gravitational 
potential for NGC 3842. Two cases are shown: constant T (a = 0; 
grey) and a T gradient with a = —0.2 (black). The black hole 
dominates the mass at small radii within the sphere of influence 
r'i n f, stars dominate at intermediate radii of r ~ 0.1r c ff, and the 
dark matter halo takes over at r > r c g . 

6086, and NGC 7768, respectively. As a increases from 
to +0.1, the best-fit values of M. increase by 14%, 
22%, and 17% for the three galaxies. For comparison, 
individual trials yield statistical errors ~ 30%-40% in 
M.. Color and line strength gradients in these three 
galaxies suggest a ~ —0.2 to —0.1. 

Our results suggest that gradients in T may be a non- 
negligible source of systematic error in measurements of 
M, and M*. The overall effect of this error on the corre- 
lations between M, and host galaxy properties will de- 
pend on the distribution of T gradients across the sample 
of galaxies with existing M m measurements. 

All three galaxies in our study are BCGs with M* > 
1O 12 M and velocity dispersions a > 250 km s -1 . Ex- 
trapolating the trend illustrated in Figure [3] to other val- 
ues of a will indicate the potential level of bias in dynam- 
ical measurements of M, that assume uniform T. A sys- 
tematic variation in a with galaxy stellar mass or velocity 
dispersion can have an impact on the slopes of the M m - 
Mbulge and M,-a scaling relations, w hich are currently 
deter mined assuming a = (e.g., iMcConnell fc Mai 
20131 and references therein). For early- type galaxies, 
Tortora et al.l (|2011l ) find a non-monotonic trend in the 
dependence of a on M* and a, where a has a mini- 
mum around M± ~ 1-3 x 10 10 M© (a ~ 100-160 kms" 1 ) 
and increases mildly toward the low- and high-M* (a) 



ends. The scatter, however, is large. If a decreases (in- 
creases) with decreasing Mbuige or a, our results here 
indicate that the inferred M, would be systematically 
lower (higher) for smaller galaxies, thereby steepening 
(flattening) the slopes of the scaling relations. In some 
galaxies where kinematic measurements superbly resolve 
ri n f, lower degeneracies between M, and T will reduce 
systematic errors in M m . 

In addition to biasing measurements of M m , gradients 
in T may impact other attempts to decompose mass pro- 
files into multiple components. In particular, strong and 
weak lensing provide complementary data to stellar kine- 
matics for probing the total mass profiles of galaxies and 
galaxy clusters. In contrast to our small-scale kinematic 
data for NGC 3842, NGC 6086, and NGC 7768, lens- 
ing studies typically constrain the total enclosed mass at 
larger radii of ~ 1 — 100 r e s- In this case, the inner slope 
of the dark matter profile can be degenerate with the 
stellar mass profile. Current methods estimate the stel- 
lar mass component using stellar kinematics and/or pop 
ulation synthe sis, under the assumption of spatially uni 
form T (e.g. Uiang fc Kochanekll2007t iSand et al.lfeoO " 
Auger et at] 120101: iTortora et al.1 l2010bt iNewman et al 
20 111 120121 ) . Allowing for a gradient in T in these stud- 



ies could affect the decomposition of the stellar and dark 
matter components. For instance, in the presence of a 
negative gradient in T, models sensitive to the enclosed 
mass slope near 1 r c g might prefer a less cuspy dark mat- 
ter profile, to counteract the inward steepening in stellar 
mass. 

Gravitational lensing is better suited for probing lens- 
ing galaxies at z ~ 0.2 and beyond, whereas direct dy- 
namical measurements of M, are currently feasible only 
in local galaxies within ~ 150 Mpc. Other dynamical 
tracers such as globular clusters have been used to in- 
fer galaxy mass profiles out to several effective radii (e.g. 
iRomanowskv et all [200llMurphv et al.l 120111) . A joint 
analysis combining the kinematics from our small-scale 
studies with those from the large-scale tracers can pro- 
vide more comprehensive measurements of the mass pro- 
files of stars and dark matter in early-type galaxies, from 
the black hole sphere of influence out to several effective 
radii. 
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